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Abstract 


Groundwater  samples  from  shallow  (<80m)  wells  in  the  Sand 
River  area  (NTS  73L)  were  collected  and  chemically  analyzed 
in  August  1974  and  in  March  1975  to  assess  the  magnitude  of 
seasonal  variation  of  groundwater  composition.  The  samples 
were  collected  along  a  west-east  traverse  that  essentially 
paralleled  the  regional  topographic  slope  and  the  inferred 
regional  groundwater  flow  direction.  The  means  and  standard 
deviations  (frequency  histograms)  of  the  concentration  dif- 
ferences (in  mg/L)  between  summer  and  winter  samplings 
were  comparable  in  magnitude  and  showed  that  no  statistical 
difference  could  be  found  between  the  two  suites  of  samples. 
The  variable  shape  and  range  of  the  frequency  distributions, 
however,  showed  that  a  simple  statistical  averaging  approach 
was  not  adequate  as  a  way  of  interpreting  the  temporal  varia- 
tions. Where  the  frequency  distributions  were  asymmetrical, 
as  in  the  case  of  sodium,  magnesium  and  chloride,  the  trend 
was  for  a  higher  summer  than  winter  concentration.  When  the 
differences  were  plotted  along  the  traverse,  generally  positive 
winter-summer  concentration  differences  were  found  in  drift 


and  bedrock  aquifers  associated  with  the  Wapiti  Formation, 
while  mainly  negative  shifts  of  the  concentration  differences 
were  found  in  the  thin  drift  aquifer  underlain  by  the  Lea  Park 
Formation.  Positive  winter  to  summer  shifts  in  the  area  of 
Wapiti  Formation  subcrop  are  likely  the  result  of  dilution  and 
mixing.  Negative  shifts,  on  the  other  hand,  were  found  primari- 
ly in  the  area  underlain  by  the  Lea  Park  Formation  and  were 
explained  as  due  to  flushing  of  salts  from  the  soil  and  un- 
saturated zones  during  summer  recharge.  The  groundwaters 
were  very  close  to  chemical  equilibrium  with  respect  to  the 
minerals  aragonite,  calcite,  dolomite,  magnesite,  siderite, 
amorphous  silica  and  Ca-montmorillonite.  These  observations 
tend  to  support  the  belief  that  chemical  variability  of  ground- 
water composition  is  constrained  by  equilibrium  reactions  with 
aquifer  minerals.  Carbonate  and  clay  minerals  appeared  to  be 
particularly  important  in  controlling  the  composition  because 
of  the  consistency  of  chemical  equilibrium  and  the  fact  that 
these  minerals  are  found  everywhere  in  the  drift  and  bedrock 
aquifers  of  the  region. 


Introduction 


The  spatial  distribution  of  cliemical  constituents  in 
groundwater  is  generally  related  to  the  flow  pattern, 
mineralogical  composition  and  climate  of  a  ground- 
water basin.  The  fact  that  certain  predictable  changes 
in  both  the  quality  and  quantity  of  dissolved  matter  take 
place  during  groundwater  movement  from  areas  of 
recharge  to  areas  of  discharge  is  useful  in  cor- 
roborating or  refining  the  physical  hydrogeological  pro- 
perties of  a  groundwater  flow  system.  In  general,  the 
chemical  composition  of  groundwater  most  closely 
resembles  the  weathering  products  of  source  materials 
in  the  recharge  area,  and  the  degree  of  mineralization 
increases  exponentially  with  depth  due  to  a  decrease 
in  flushing  intensity,  and  laterally  due  to  further 
dissolution  of  minerals  along  the  flow  path. 

Temporal  variations  of  groundwater  chemistry 
reflect  hydrogeochemical  processes  operating  near 
the  sampling  point.  Such  processes  can  be  both 
physico-chemical  and  biochemical,  and  the  magnitude 
of  variation  in  concentration  of  chemical  constituents 
may  be  great  enough  to  contribute  uncertainty  to  the 
variability  attributed  to  the  spatial  distribution.  In  this 
paper,  the  chemical  variation  between  groundwater 
samples  collected  during  winter  and  summer  was  a 
reasonable  first  attempt  at  documenting  the 
phenomenon  of  temporal  variation  on  a  recon- 
naissance scale. 1 

This  paper  documents  spatial  and  temporal  varia- 
tions in  groundwater  composition  and  shows  how 
these  variations  relate  to  regional  hydrogeological 
characteristics  on   a  reconnaissance  scale  of 

1  Davison  and  Vonhof  (1978)  have  reported  on  groundwater  chemical 
variations  within  an  area  of  approximately  12  km2  associated  with  a 
buried  channel  aquifer  near  Esterhazy,  Saskatchewan.  Samples 
were  taken  during  November  and  June,  July  and  August  at  intervals 
of  one  day. 


1 :250  000.  The  role  of  solution-mineral  equilibrium  as  a 
stabilizing  influence  on  groundwater  chemical  com- 
position is  also  considered. 

The  study  area  encompasses  some  14  000  km2  and 
is  about  95  km  northeast  of  Edmonton  (see  figure  1). 
The  area  lies  entirely  within  the  forested  Interior  Plains 
physiographic  province.  Highest  elevations  (about  860 
m  amsl)  are  found  in  the  north-central  part  of  the  area; 
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Figure  1.  Location  map. 
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the  Beaver  River  valley  at  the  point  of  juncture  with  the  mean  annual  precipitation  ranges  from  363  to  620  mm. 
southeastern  boundary  of  the  area  has  an  elevation  of  The  area  is  almost  entirely  within  the  "short,  cool  sum- 
490  m.  Mean  annual  temperature  is  about  1  °C,  and  the     mer"  climatic  zone  (Ozoray  et  al.,  1980). 


Physical  hydrogeology 


The  most  productive  bedrock  aquifer  in  the  study  area 
is  the  upper  Cretaceous  Wapiti  Formation,  but  the 
lithologies  alternate  between  lenticular  mont- 
morillonitic  sandstones  and  clay  shales,  and  the 
groundwater  yields  vary  between  25  to  100  L/min. 
Generally,  the  glacial  drift  (consisting  of  sand,  sand 
and  silt  or  gravel,  and  till)  is  the  most  common  shallow 
aquifer,  and  is  characterized  in  most  detail  in  this 
report.  However,  the  variability  of  drift  lithology  and 
thickness  (as  little  as  3  m  in  the  southwest  highland  to 
more  than  130  m  in  parts  of  buried  valleys)  is  reflected 
in  a  range  of  yields  from  5  to  25  L/min  in  some  silt  or  till 


deposits  to  450  to  2250  L/min  in  some  sand  and  gravel 
lenses.  The  Lea  Park  Formation,  also  upper 
Cretaceous  in  age,  is  a  poorly  permeable  marine  clay- 
shale;  it  is  the  major  bedrock  aquitard  in  the  area  and 
the  top  of  this  formation  essentially  marks  the  lower 
limit  of  groundwater  supplies  in  the  area.  Groundwater 
availability  is  very  poor  where  the  drift  cover  over  the 
Lea  Park  is  thin.  The  regional  flow  of  groundwater  is 
from  west  to  east,  and  is  essentially  parallel  to  the 
trends  of  major  streams  and  buried  channels  (Ozoray 
etal.,  1980). 


Chemical  hydrogeology 


Acquisition  of  existing  chemical  data 

Chemical  analyses  of  groundwater  sampled  from  wells 
in  drift  and  shallow  bedrock  aquifers  were  obtained 
from  the  Central  Data  File  of  the  Groundwater  Depart- 
ment, Alberta  Research  Council.  The  analyses  were, 
for  the  most  part,  the  work  of  the  Provincial  Analyst  and 
the  Pollution  Control  Laboratory  of  Alberta  Environ- 
ment. The  analyses  were  tested  for  internal  precision 
by  means  of  a  cation-anion  balance  (where  a  complete 
analysis  was  given  of  major  constituents)  or  a  total 
dissolved  solids  (TDS)  balance  when  determinations  of 
sodium,  potassium  and  separate  analyses  of  calcium 
and  magnesium  were  lacking.  Of  the  588  groundwater 
chemical  analyses  that  were  coded,  183  mainly  partial 
analyses  were  rejected  because  error  (e)  exceeded  20 
percent  for  the  TDS  balance.  The  expression  that  was 
evaluated  is  written  below: 


e  (o/o)  = 


cations 


-HE. 


TDS 


TDS 


X  100 


*Eanions  =  SO4  +  C\  +  NO3  +  CO3  +  0.49  HCO3 

Field  sampling  and  analysis 

An  additional  90  groundwater  samples  were  first  col- 
lected along  Highway  28  from  Vilna  to  Cold  Lake  in 
August  1974  and  a  second  sampling  of  the  same  wells 
took  place  in  March  1975.  The  path  of  the  traverse  and 
the  distribution  of  wells  sampled  are  shown  on  figure  2. 
The  path  cut  across  the  area  of  greatest  data  density 
and  was  generally  parallel  to  the  regional  west  to  east 
topographic  slope  in  the  southern  half  of  the  map  area: 
from  a  high  in  the  southwest  of  670  masi  near  Vilna  to  a 


low  at  Cold  Lake  of  533  masl.  Data  of  sufficient  density 
for  study  were  found  only  in  the  southern  half  of  the 
map  area  and  in  the  Lac  La  Biche  Reserve.  The 
traverse  is  also  assumed  to  be  parallel  to  any  regional 
flow  components. 

The  groundwater  samples  were  collected  from  pro- 
ducing domestic  wells  after  pumping  until  constant  pH 
and  temperature  readings  were  obtained  (usually 
within  15  minutes).  Three  2-litre  samples  were  col- 
lected in  polyethylene  bottles:  one  for  field  measure- 
ment of  electrical  conductivity,  pH,  alkalinity,  iron  con- 
tent and  temperature;  a  second  acidized  (drop  by  drop 
addition  of  concentration  HNO3  to  pH  <2.0)  for 
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analysis  of  total  iron,  calcium  and  magnesium  con- 
tents; and  a  third  for  determination  of  major  consti- 
tuents. The  latter  two  samples  were  analyzed  using 
standard  methods  in  the  Geological  Survey  Chemistry 
Laboratory  of  the  Alberta  Research  Council. 

Measurement  of  field  pH  was  carried  out  using  an 
ORION  model  401  meter  and  ORION  91-92  combina- 
tion electrode.  Buffers  (pH  =  7.00  and  pH  =  4.00) 
were  refrigerated  at  about  10°C,  a  temperature  ap- 
proximating that  of  the  groundwater.  Alkalinity  was 


determined  by  titrating  50  mL  of  water  with  0.01647  N 
HCI  to  obtain  pH  versus  volume  of  titrant  added.  The 
method  of  Barnes  (1964)  was  employed  to  calculate 
the  field  alkalinity.  Conductivity  was  measured  using  a 
Hach  model  2510-01  conductivity  meter.  A  Hach  model 
2504  spectrophotometer  was  used  to  analyze  the  total 
iron  content  by  means  of  the  1,10  -  Phenanthroline 
method.  Temperature  of  the  groundwater  was 
measured  with  a  glass  mercury  thermometer 
calibrated  in  0.2°C  intervals. 


Results 


General  aspects  of 
groundwater  chemistry 

High  TDS  (<1000  mg/L)  are  found  in  groundwater 
above  or  within  the  Wapiti  Formation;  groundwater 
from  buried  channel  aquifers  that  overlie  the  Lea  Park 
marine  shale  generally  tends  to  be  less  than 
1000  mg/L.  Figure  3  shows  these  observations.  Note 
that  local  pockets  of  either  very  poor  (>2000  mg/L)  or 
very  fresh  (<500  mg/L)  groundwater  occurs  within  the 
area. 

As  shown  in  figure  4,  the  most  abundant  chemical 
type  of  groundwater  is  a  mixed  type  with  neither  bicar- 
bonate nor  sulfate  type  being  dominant. ^  Regions  of 
calcium-magnesium-bicarbonate  groundwater  mainly 
occur  in  the  more  productive  drift  aquifers.  Calcium- 
sulfate  and  sodium-sulfate  type  groundwaters  occur 
primarily  in  the  southwestern  portion  of  the  area. 

Groundwaters  from  the  area  are  generally  very  low  in 
chloride,  although  concentrations  as  high  as  100  mg/L 
were  found  in  the  northwestern  corner  of  the  Cold  Lake 


Indian  Reserve.  Near  Goodridge,  concentrations  as 
high  as  70  mg/L  were  noted.  Chloride  concentrations 
are  generally  well  below  50  mg/L  elsewhere  in  the 
Sand  River  area,  and  are  especially  low  in  aquifers 
associated  with  the  Wapiti  Formation  and  the  pre- 
glacial  channels,  presumably  due  to  active  circulation 
of  groundwater. 

The  distribution  of  groundwater  chemical  types  with 
respect  to  depth  can  be  portrayed  by  a  modified  Piper 
diagram  (Piper,  1944).  The  mean  depth  of  the  open  in- 
terval is  plotted  in  the  water-type  field  in  figure  5.  The 
figure  shows  that  calcium-magnesium-bicarbonate 
waters  predominate  at  depths  less  than  30  m,  although 
other  water  types  still  appear  in  this  depth  zone. 
Calcium-magnesium-sulfate  type  groundwater  is 
generally  found  at  depths  less  than  45  m.  Sodium- 
sulfate  type  of  groundwater  tends  to  occur  in  deeper 
wells.  The  trend  indicates  the  groundwaters  change 
from  calcium-magnesium-bicarbonate  facies  to 
sodium-sulfate  facies  with  increasing  depth. 


4 


0         10       20       30       40       50       60      70       80  90 
S04=  +  C|-  +  N03-(epm%) 

Figure  5.  Distribution  of  groundwater  chemical  types  with 
respect  to  depth. 


Results  of  field  study 

Table  1  summarizes  the  results  of  the  field  and 
laboratory  chemical  analyses.  Chemical  analyses  of 
major  and  minor  constituents  were  performed  in  the 
Geological  Survey  Chemistry  Laboratory,  under  the 
supervision  of  J.R.  Nelson,  within  two  or  three  weeks 
after  sampling.  The  symbols,  Ca^**  and  Mg^^*  are  the 
concentrations  of  calcium  and  magnesium  respective- 
ly in  the  acidified  samples.  The  adjusted  bicarbonate 
concentration,  HCOsy,  takes  the  precipitation  of  calcite 
between  times  of  sampling  and  analysis  into  account 
(Wallick,  1977).  The  table  gives  the  results  of  summer 
and  winter  sampling. 

Some  observations  concerning  the  magnitude  and 
variation  of  the  chemical  parameters  may  be  made: 

Differences  between  the  concentrations  of  various 
constituents  from  winter  and  summer  samplings  are 
expressed  as  frequency  histograms  (figure  6  a-h).  In  all 


cases,  the  mean  is  not  significantly  different  from  zero, 
indicating  no  clear  difference  between  winter  and  sum- 
mer groundwater  composition.  However,  the  variable 
shapes  and  ranges  of  the  frequency  histograms  reflect 
an  aspect  of  the  chemical  variability  that  is  not  explain- 
ed by  a  simple  statistical  averaging  approach.  For  ex- 
ample, A[Ca^*]  is  distributed  more  or  less  symmetrical- 
ly about  zero  as  a  bell-shaped  normal  distribution. 
A  [Na*],  however,  is  not  symmetrically  distributed,  but 
skewed  toward  negative  differences.  Examination  of 
the  remaining  distributions  shows  that  AfHCOa"], 
A[S04=],  A[N031,  and  ATDS  are  normal,  while  A[Mg**] 
and  A[CI"]  are  skewed.  The  sodium,  magnesium  and 
chloride  distributions  are  skewed  in  the  direction  of 
higher  summer  than  winter  concentrations. 

Another  way  of  looking  at  the  differences  is  to  con- 
sider their  variation  with  distance  along  the  traverse. 
Figure  7  (a-h)  shows  the  differences  expressed  as  a 
percent  of  the  summer  concentration.  The  diagrams 
demonstrate  the  relative  importance  of  the  variations, 
showing  that  in  excess  of  100  percent  variation  in  Na* 
can  be  obtained  and  more  than  50  percent  variation  in 
Ca**  and  Mg*^  occurred.  With  respect  to  anions,  the 
variation  in  HCO3"  was  usually  less  than  50  percent, 
[CI"]  showed  large  percentage  variation  (in  excess  of 
100  percent  for  positive  shifts  and  on  the  order  of  50 
percent  for  negative  shifts),  and  [S04=]  and  [NO3"] 
showed  similarly  large  variations. 

Plotting  these  variations  along  the  traverse  is  a  sim- 
ple means  of  checking  whether  the  differences  are  ran- 
domly distributed  in  space.  They  are  not.  Positive 
deflections  of  Na*,  Mg**  and  S04=  tend  to  occur  along 
the  west  half  of  the  traverse  while  negative  deflections 
tend  to  occur  along  the  east  half.  Variations  of  CI"  and 
NO3"  generally  tend  to  be  uniform  along  the  traverse. 


Table  1.  Chemical  analyses  of  groundwater  samples,  Sand  River  (73L) 


Location 

Date 

Depth 

Temp 

PH, 

PHl 

Cond 

Ca*^ 

Ca-T 

Na* 

CO3  HCO3-T 

4-61-6-14-1 

24/8/74 

4.3 

8.89 

7.50 

8.20 

1150 

73.0 

108.7 

73.5 

82.8 

16.3 

3.3 

-  512.6 

4-61-6-14-1 

25/3/75 

4.3 

9.00 

7,10 

7.50 

990 

102.0 

113.3 

67.0 

71.1 

18.8 

3.3 

-  654.0 

4-61-5-22-13 

25/8/74 

5.5 

12.22 

7.32 

8.00 

2700 

116.0 

142.9 

300.0 

300.0 

78.8 

5.0 

-  1022.7 

4-61-5-22-13 

24/3/75 

5.5 

3.50 

7.30 

7.60 

1610 

137.0 

166.7 

137.0 

133.3 

41.3 

9.6 

-  576.0 

4-61-5-26-12 

25/8/74 

6.7 

1.22 

7.12 

8.00 

2000 

119.2 

156.2 

195.0 

194.4 

89.0 

5.8 

-  905.5 

4-61-5-26-12 

25/3/75 

6.7 

7.00 

7.20 

7.40 

1025 

100.0 

114.4 

63.0 

67.8 

27.5 

5.8 

-  625.0 

4-61-5-25-1 

25/8/74 

30.5 

8.33 

7.35 

7.90 

4000 

279.0 

322.7 

325.0 

350.0 

399.0 

8.3 

-  571.1 

4-61-5-25-1 

24/3/75 

30.5 

8.80 

7.30 

7.80 

2600 

187.0 

221.1 

78.0 

80.0 

289.0 

8.8 

-  598.0 

4-61-4-32-1 

25/8/74 

9.8 

8.89 

7.32 

8.20 

1000 

77.7 

77.0 

66.0 

64.4 

64.0 

3.3 

-  617.5 

4-61-4-32-1 

24/3/75 

9.8 

4.00 

7.30 

7.60 

980 

76.0 

88.9 

63.0 

67.8 

55.0 

4.6 

-  615.0 

4-62-4-2-8 

25/8/74 

8.5 

9.44 

7.30 

7.60 

850 

87.6 

91.3 

56.0 

57.2 

11.0 

2.9 

-  524.8 

4-62-4-2-8 

24/3/75 

8.5 

7.30 

7.50 

940 

96.0 

108.9 

64.0 

66.7 

11.3 

4.2 

-  586.0 

4-63-2-24-9 

26/8/74 

17.1 

11.66 

7.03 

7.80 

1700 

172.6 

186.2 

99.5 

95.6 

109.0 

6.3 

-  588.2 

4-63-2-24-9 

24/3/75 

17.1 

11.00 

6.90 

7.30 

1880 

199.0 

225.6 

99.0 

104.4 

105.0 

9.2 

-  615.0 
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Table  1.  Chemical  analyses  of  groundwater  samples,  Sand  River  (73L)  (cont'd) 


Location 

Date 

Depth 

Temp 

pHf 

PHl 

Cond 

Ca** 

Ca**-p 

Na* 

CO3  HCO3 

4-63-2-15-1 

IQIQIIA 

12.2 

8.33 

7.50 

7.90 

1575 

130.0 

127.8 

108.5 

103.9 

81 .0 

3.8 

624.8 

4-63-2-15-1 

24/3/75 

12.2 

3.80 

7.10 

7.70 

965 

1 10.0 

1 24.4 

56.0 

58.9 

25.0 

5.8 

-  493.0 

4-62-2-35-15 

26/3/74 

1 1 .4 

8.89 

7.07 

8.00 

1000 

101 .8 

1 13.1 

51 .0 

47.2 

34.0 

6.7 

-  449.1 

4-62-2-35-15 

24/3/75 

1 1 .4 

5.50 

6.90 

7.30 

815 

105.0 

123.3 

48.0 

53.3 

23.8 

7.9 

-  454.0 

4-62-2-36-14 

26/3/74 

9.1 

14.44 

7.1 8 

8.00 

1 100 

154.8 

1 52.0 

52.0 

52.2 

14.0 

4.2 

546.7 

4-62-2-36-14 

24/3/75 

9.1 

5.50 

7.00 

7.40 

825 

1 06.0 

1 1 7.8 

47.5 

50.0 

7.5 

1 .7 

527.0 

4-62-1-31-15 

26/8/74 

7.3 

8.89 

7.23 

8.30 

8.50 

92.8 

96.4 

48.5 

48.9 

1 1 .0 

5.4 

-  468.6 

4-62-1-31-15 

24/3/75 

7.3 

7.30 

7.10 

7.50 

870 

101 .0 

111.1 

47.5 

53.3 

6.3 

5.0 

503.0 

4-63-1-2-2 

26/8/74 

7.6 

6.1 1 

7.20 

7.40 

950 

122.2 

129.6 

50.0 

50.0 

5.0 

2.9 

566.3 

4-63-1-2-2 

24/3/75 

7.6 

9.00 

7.40 

8.10 

875 

88.0 

1 24.4 

45.6 

46.7 

5.0 

1 0.4 

464.0 

4-62-2-22-16 

27/8/74 

29.0 

9.44 

7.20 

7.80 

1 100 

133.2 

136.0 

57.5 

55.0 

54.0 

5.0 

663.9 

4-62-2-22-16 

24/3/75 

29.0 

4.70 

71 .0 

8.20 

1 1 20 

111.0 

1 36.7 

55.0 

56.7 

40.0 

2.9 

605.0 

4-62-2-23-4 

27/8/74 

1 1 .6 

8.89 

6.97 

7.60 

1 100 

1 59.8 

1 66.7 

55.5 

57.8 

15.0 

63. 

541.9 

4-62-2-23-4 

24/3/75 

1 1 .6 

6.00 

7.00 

7.20 

1110 

1 74.0 

1 93.3 

53.0 

55.6 

7.5 

5.8 

542.0 

4-62-2-14-5 

27/8/74 

10.4 

14.44 

6.91 

7.90 

875 

1 13.4 

129.1 

38.5 

40.0 

1 5.0 

7.0 

-  488.2 

4-62-2-14-5 

10.4 

1 1 .0 

6.90 

7.30 

910 

134.0 

145.6 

39.7 

42.2 

8.8 

6.3 

517.0 

4-62-3-8-13 

27/5/74 

5.5 

6.1 1 

7.07 

7.60 

3300 

253.3 

261 .7 

1 78.0 

1 66.7 

356.0 

4.6 

868.9 

4-62-3-8-13 

24/3/75 

5.5 

4.50 

7.10 

7.80 

2500 

148.0 

208.9 

1 07.0 

1 07.8 

258.0 

5.4 

615.0 

4-59-13-20-1 

21/8/74 

15.2 

3.89 

7.40 

8.10 

1600 

63.7 

1 58.3 

68.0 

69.9 

68.8 

22.9 

361.2 

4-59-13-20-1 

OC  /o  /^c 

26/3/75 

15.2 

1 .30 

7.30 

7.60 

1050 

1 19.0 

1 36.7 

48.5 

50.0 

52.5 

6.3 

608.0 

4-59-1  d-iiO-b 

21/8/74 

39.6 

2.22 

7.40 

8.10 

1400 

98.4 

1 29.6 

53.0 

54.4 

55.0 

4.6 

478.4 

A    CC\    H  O   OO  C 

4-59-13-22-5 

OC  /O  /TC 

26/3/75 

39.6 

4.50 

7.30 

7.50 

1210 

138.0 

1 58.9 

54.0 

57.8 

61 .0 

5.0 

603.0 

y1   CO  H  O  OO  C 

4-59-1  o-2ii-5 

21/8/74 

42.7 

4.44 

7.48 

7.80 

900 

89.5 

96.0 

38.5 

41 .1 

1 8.8 

6.7 

466.2 

y1    CO   H  O   OO  C 

4-59-13-33-5 

OC  /o  /^c 

26/3/75 

42.7 

3.00 

7.30 

7.60 

790 

93.0 

106.7 

35.0 

38.9 

23.8 

6.7 

-  454.0 

4-59-13-33-2 

22/8/74 

i  5.2 

13.88 

7.02 

7.70 

1 100 

155.4 

1 61 .1 

61 .5 

64.4 

8.8 

3.3 

-  732.2 

A   CO  H  O  OO  O 

4-59-12-32-2 

OO  /C  /"7C 

23/6/75 

5.2 

6.80 

7.00 

7.30 

1300 

180.0 

203.3 

62.0 

63.3 

1 5.0 

7.1 

571 .0 

4-59-12-33-1 

22/5/74 

30.5 

9.44 

7.20 

7.60 

1450 

170.3 

180.9 

64.5 

64.9 

98.8 

6.3 

-  639.5 

A    CO   H  O   OO  H 

4-59-12-33-1 

25/3/75 

30.5 

5.70 

7.20 

7.60 

1550 

173.0 

196.7 

66.0 

68.9 

104.0 

6.7 

644.0 

4-59-1 1-28-12 

OO  /O  Z"?  jl 

22/8/74 

24.4 

6.10 

7.29 

7.70 

1550 

156.0 

158.1 

106.0 

111.1 

16.5 

6.7 

-  656.6 

>l    CO    H  ^    OO    H  O 

4-59-1 1-28-12 

25/3/75 

24.4 

2.60 

7.10 

7.40 

1820 

214.0 

243.3 

103.0 

1 06.7 

46.3 

6.7 

-  561.0 

A    CC\    H  ^    OO   H  ^ 

4-59-1 1-30-16 

OO  /O  >l 

22/8/74 

15.2 

6.70 

7.35 

7.80 

1000 

1 10.8 

1 15.1 

54.5 

57.8 

1 .3 

5.8 

-  556.5 

A   CO   H  H    OO   ^  C 

4-59-1 1-30- lb 

25/3/75 

15.2 

4.50 

7.20 

7.50 

990 

144.0 

215.6 

59.0 

76.7 

12.5 

10.0 

-  564.0 

>l    CO    ^  H    OO    H  ^ 

4-59-1 1-28-16 

OO  /O  >! 

22/8/74 

35.4 

7.20 

7.22 

7.60 

1850 

185.0 

203.1 

78.0 

82.2 

10.13 

5.4 

-  676.1 

>1   CO   -4  H    OO   ^  C 

4-59-1 1-28-16 

25/3/75 

35.4 

2.50 

7.30 

7.50 

1700 

182.0 

133.3 

71 .0 

60.0 

1 13.0 

5.4 

639.0 

A  cr\  H  o  o  o 

4-60-10-8-8 

22/8/74 

8.5 

6.10 

7.18 

7.50 

1700 

197.0 

221 .7 

108.0 

108.9 

28.8 

6.3 

-  693.0 

A   CO  A  f\  Q  A 

4-60-10-8-4 

25/3/75 

8.5 

3.40 

7.20 

7.70 

1925 

220.0 

253.3 

124.0 

1 15.6 

36.3 

6.7 

800.0 

A    Cf\   A  O  OC  H 

4-60-10-25-1 

OO  /O  i1 

23/8/74 

76.2 

7.80 

7.57 

8.10 

2150 

72.6 

88.9 

28.5 

30.0 

431 .0 

5.0 

-  551.6 

A   CO  H  O  OC  H 

4-60-10-25-1 

25/3/75 

76.2 

5.00 

7.50 

7.70 

2500 

93.0 

104.4 

31 .0 

32.2 

426.0 

5.4 

-  600.0 

>4   CO  O  H  *7  O 

24/8/74 

12.2 

7.20 

7.46 

8.20 

600 

63.8 

69.2 

25.0 

26.7 

5.0 

2.5 

-  346.6 

A  c^r\  O  H  "7  o 

4-b(J-y- 1  f-d 

25/3/75 

12.2 

5.00 

7.30 

7.70 

625 

71 .0 

78.9 

25.0 

26.7 

283.0 

3.8 

371.0 

A   CO  O  H  O  O 

24/8/74 

59.7 

7.20 

7.80 

8.20 

2075 

39.2 

42.3 

16.0 

15.6 

427.5 

2.9 

-  634.6 

A  CO  Q  H  O  H  O 

25/3/75 

59.7 

3.80 

7.70 

7.60 

2020 

42.0 

56.7 

16.8 

16.7 

444.0 

4.6 

-  637.0 

A   CO  C   H  C  >1 

4-60-5-16-4 

24/5/74 

16.5 

6.10 

7.1 1 

7.60 

1700 

145.0 

143.3 

128.0 

137.8 

18.8 

4.2 

-  595.6 

A    Cf\   C    A  C  A 

4-60-5-16-4 

25/3/75 

16.5 

4.50 

7.20 

7.20 

1880 

164.0 

162.0 

139.0 

147.0 

27.5 

4.6 

-  639.0 

4-60-7-1 7-1 

24/5/74 

12.2 

8.90 

7.10 

7.90 

1650 

139.0 

178.9 

105.0 

108.3 

30.0 

5.0 

-  532.1 

4-60-7-17-1 

25/3/75 

12.2 

7.50 

7.10 

7.30 

1550 

161 .0 

168.9 

92.0 

102.2 

37.5 

5.4 

-  610.0 

4-60-7-12-15 

24/8/74 

10.0 

6.10 

7.25 

7.80 

1100 

120.0 

128.8 

48.0 

51.1 

28.8 

5.0 

-  649.3 

4-60-7-12-16 

25/3/75 

10.1 

3.80 

7.10 

7.30 

1040 

120.0 

144.4 

51.0 

55.6 

36.3 

5.0 

-  676.0 

4-60-6-16-7 

24/8/74 

30.0 

8.90 

7.21 

7.80 

1700 

102.7 

150.6 

110.0 

124.4 

56.3 

5.8 

-  561.4 

4-60-6-16-9 

28/5/74 

30.0 

4.80 

7.00 

7.30 

1210 

163.0 

185.6 

64.0 

65.6 

6.3 

9.6 

-  698.0 

4-60-6-35-13 

24/5/74 

5.5 

5.60 

7.15 

7.90 

1200 

107.0 

136.9 

47.5 

52.8 

22.5 

6.7 

-  451.5 

4-60-6-35-13 

25/3/75 

5.5 

2.30 

7.20 

7.60 

925 

96.0 

115.6 

48.0 

53.3 

48.8 

4.2 

-  617.0 
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Figures  6  a-h.  Concentration  differences  between  winter  and  summer  sampling  of  the  same  wells,  expressed  in  frequency 
histograms. 
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Figures  7a-h.  Concentration  differences  between  winter  and  summer  sampling  of  the  same  wells  along  the  traverse,  express- 
ed as  mg/L  percent. 


Discussion 


The  point  where  the  traverse  crosses  the  contact  bet- 
ween the  Wapiti  Formation  and  the  Lea  Park  Forma- 
tion is  indicated  on  figures  7  (a-h)  by  a  vertical  line.  It 
may  be  noted  that  this  point  correlates  with  the 
changes  noted  in  the  previous  section  pertaining  to 
positive  versus  negative  percentage  differences.  On 


figure  3,  a  contour  map  of  TDS  is  superimposed  upon  a 
map  giving  the  contact  between  the  Belly  River  and 
Lea  Park  Formations  and  the  bedrock  topography. 
Note  that  the  TDS  are  significantly  higher  in  ground- 
water obtained  from  aquifers  in  the  area  of  Wapiti  For- 
mation subcrop  compared  to  the  TDS  from  aquifers 
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overlying  the  Lea  Park  Formation.  This  difference  oc- 
curs for  the  following  reasons:  the  aquifer  thickness 
and  permeability  are  greater  in  the  area  underlain  by 
the  Wapiti  Formation  than  in  the  area  underlain  by  the 
Lea  Park,  where  the  aquifers  are  thin  drift  sands  and 
channel  sands  and  gravels  resting  on  very  poorly 
permeable  marine  clay-shale.  Therefore,  deeper  flow 
systems  with  longer  groundwater  residence  times  and 
greater  groundwater  mineralization  can  develop  in  the 
area  associated  with  the  Wapiti  Formation. 

Positive  winter-summer  concentration  differences  in 
aquifers  associated  with  the  drift  and  bedrock  in  the 
area  of  Wapiti  Formation  subcrop  may  also  be  ex- 
pected because  of  the  effects  of  dilution  and  mixing  as 
a  result  of  groundwater  recharge.  On  the  other  hand, 
negative  fluctuations  of  the  percentage  concentration 
differences  in  the  thin  drift  aquifer  underlain  by  the  Lea 
Park  ( see  particularly  figure  7 )  can  be  explained  by 
flushing  of  salts  from  the  soil  and  unsaturated  zones 
during  summer  recharge.  The  total  volume  of  ground- 
water in  storage  in  the  thin  drift  aquifer,  which  is 
significantly  less  than  the  volume  stored  in  the  thicker 
drift  and  bedrock  aquifer  associated  with  the  Wapiti 
Formation,  suggests  that  the  effect  of  flushing  of  salts 
would  be  accentuated  in  aquifers  overlying  the  Lea 
Park.  The  relative  effect  of  any  dilution  or  flushing 
event  would  appear  greater  in  aquifers  that  are  of  small 
storage  capacity  because  of  the  greater  proportion  of 
added  salt  or  water  to  that  already  in  storage. 

Control  of  solubility  by  minerals 

Whereas  the  time  and  spatial  variability  of  the  chemical 
concentrations  in  groundwater  are  increased  by  mixing 
and  flushing,  equilibrium  controls  on  solubility  by 
aquifer  minerals  tend  to  stablize  the  concentrations.  It 
is  useful,  therefore,  to  determine  the  disequilibrium  in- 
dex distribution  with  respect  to  various  minerals  pre- 
sent in  the  drift  and  bedrock  of  the  Sand  River  area. 

Relation  of  hydrochemistry 
to  aquifer  mineralogy 

The  disequilibrium  index  was  defined  by  Paces  (1972) 
as  the  log  of  the  ratio  of  the  activity  product  to  the 
temperature-corrected  equilibrium  constant: 
I  =  log  (Q/K^) 

For  example,  in  the  case  of  gypsum  dissolved  in  pure 
water,  the  chemical  reaction  would  be 

CaS04*2H20  ^  Ca"*  +  864=  +  2H2O 

The  activity  product  for  the  above  reaction  is 

(Ca-)(S04=) 

(CaS04) 

where  ( )  =  free  ion  activity 
If  the  solution  were  undersaturated  or  oversatu rated, 
the  value  of  I  would  be  less  than  or  greater  than 


±  0.5  respectively  (Paces,  1972).  Water  samples  that 
are  saturated,  or  in  equilibrium  with  a  given  mineral, 
are  assumed  to  have  been  in  contact  with  that  mineral 
prior  to  pumping  or  spring  discharge. 

Disequilibrium  indices  with  respect  to  calcite, 
dolomite,  gypsum,  Ca-montmorillonite,  illite,  kaolinite, 
albite,  anorthite  and  microcline  were  computed  by 
means  of  a  simple  computer  program.  The  program 
calculates  the  equilibrium-free  ion  activities  of  Ca**, 
Mg*^  Na^  K*,  H^  S04=,  H2CO3,  HCO3-,  C03=,  Cl-  and 
PcQ^(atm.)  by  iteratively  forming  the  ion  pairs:  CaS04, 
NaS04",  MgS04  and  CaCOs.  Input  to  the  program  con- 
sists of  the  field  pH  and  temperature,  and  bicar- 
bonate/carbonate concentration(s)  measured  in  the 
field  or  reconstructed  from  suitable  laboratory  data 
(Wallick,  1977),  Ca**  and  Mg**  measured  on  an 
acidified  portion  of  the  groundwater  sample,  and 
analyses  for  S04=,  CL",  NO3",  H4Si04,  F",  and  total  iron 
measured  on  the  acidified  sample  or  the  equivalent 
analysis  in  the  field.  The  program  was  used  to  compute 
disequilibrium  indices  for  samples  collected  in  1974. 

SOLMNEQ  (Kharaka  and  Barnes,  1973),  the  solution 
mineral  equilibrium  program,  was  used  to  compute 
disequilibrium  indices  for  laboratory  and  field  chemical 
states  of  all  samples  taken  in  1974  and  1975.  Values  of 
I  for  the  minerals  aragonite,  calcite,  dolomite,  gypsum, 
magnesite,  siderite  and  amorphous  silica  were  closest 
to  saturation  from  among  157  mineral  species  treated 
in  the  program. 

Data  from  both  computer  calculations  are  summariz- 
ed in  table  2.  Assuming  the  -0.5  <  I  <  -i-  0.5  criterion 
to  denote  equilibrium  conditions,  the  groundwater  was, 
on  the  average,  in  equilibrium  with  calcite,  aragonite, 
dolomite,  magnesite,  siderite,  amorphous  silica  and 
Ca-montmorillonite.  Of  interest  are  the  differences  in 
index  values  between  field  and  laboratory  chemical 
states  for  the  carbonate  minerals.  (The  laboratory 
chemical  state  is  the  chemical  analysis  of  the  non- 
acidified  sample  at  ~25°C  and  includes  the  laboratory 
pH.)  In  all  cases  (except  siderite),  there  is  a  shift  from  a 
saturated  to  an  oversaturated  state.  The  primary 
reason  for  this  is  the  loss  of  dissolved  CO2  from  warm- 
ing the  sample.  For  instance,  the  average  partial 
pressures  of  CO2,  Pco2'  calculated  using  SOLMNEQ, 
were  as  follows: 

a.  Pco2  (field  samples  August/74)  = 
(3.13  ±  1.89)  X  10-2  atm 

b.  Pco^  (field  samples  March/75)  = 
(2.97  ±  1.29)  X  10-2  atm 

c.  Pco,  (field  all  samples)  =  (3.08  ±  1.22)  x  10-2atm 

d.  Pco' (lab  all  samples)  =  (1.28  ±  .52)  x  10-2atm 

These  data  show  that  the  equilibrium  partial 
pressure  of  CO2  is  constant  under  field  conditions  from 
summer  to  winter,  but  that  more  than  50  percent  of  the 
CO2  escapes  by  the  time  the  sample  is  analyzed  in  the 
laboratory.  Table  2  further  documents  the  effect  of  CO2 
loss  by  sample  warming.  Note  how  close  to  saturation 
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(0  within  ±  0.5)  the  groundwaters  were  in  the  field  with 
respect  to  the  carbonate  minerals:  calcite,  aragonite, 
dolomite,  siderite  and  magnesite;  and  how,  in  the 
laboratory  state,  the  waters  were  oversatu rated.  In 
general,  the  Sand  River  area  groundwaters  analyzed 
for  this  report  were  in  equilibrium  with  carbonate 
minerals  in  the  field.  The  values  for  gypsum  in  the  field 
and  in  the  laboratory  are  identical  because  no  loss  of 
S04=  occurs  during  sample  warming. 

With  regard  to  the  silicate  minerals,  the  ground- 
Table  2.  Disequilibrium  indices,  Sand  River  (73L) 


1 

N 

Program 

Sta 

Calcite  (CaCOa) 

0.25 

0.16 

54 

W 

F 

0.17 

0.13 

88 

S 

F 

0.89 

0.17 

88 

s 

L 

Aragonite  (CaCOa) 

0.14 

0.13 

88 

s 

L 

0.78 

0.17 

88 

s 

L 

Dolomite  (CaMg(C03)2) 

0.51 

0.37 

54 

w 

F 

0.37 

0.52 

87 

s 

F 

1.70 

0.67 

87 

s 

L 

Siderite  (FeCOs) 

0.12 

0.38 

56 

s 

F 

No  laboratory  calculation  possible  because  iron 
precipitates  from  solution. 


Gypsum  (CaS04  2H2O) 

-1.46 

0.49 

54 

W 

F 

-1.21 

0.41 

89 

S 

F 

-1.22 

0.39 

87 

S 

L 

Magnesite  (MgCOa) 

-0.21 

0.20 

89 

S 

F 

0.46 

0.21 

89 

S 

L 

Amorphous  silica 

-0.12 

0.06 

88 

S 

F 

-0.24 

0.06 

85 

S 

L 

Ca-montmorillonite 

0.92 

0.10 

54 

W 

F 

lllite 

-0.73 

0.20 

54 

W 

F 

Kaolinite 

0.92 

0.16 

54 

W 

F 

Microcline 

-1.49 

0.27 

54 

W 

F 

Albite 

-2.74 

0.60 

54 

W 

F 

Anorthite 

-10.79 

0.31 

54 

W 

F 

N  =  Number  of  samples 

F  =  Field 

L  =  Laboratory 

W  =  Equilibrium  computer  model  PHCALC 
S  =  SOLMNEQ 


waters  were  in  equilibrium  with  amorphous  silica  and 
Ca-montmorillonite,  slightly  oversaturated  with  respect 
to  kaolinite,  undersaturated  with  respect  to  illite, 
microcline  and  albite,  and  greatly  undersaturated  with 
respect  to  anorthite. 

Another  way  of  interpreting  the  results  pertaining  to 
the  silicate  minerals  in  table  2  is  the  aspect  of  mineral 
stability  in  the  aquifer.  The  more  negative  the  value  for 
I,  the  less  stable  a  particular  mineral  is  when  in  contact 
with  the  specified  groundwater  composition.  Anorthite 
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is,  therefore,  the  least  stable  mineral  in  the  table,  and 
Ca-montmorillonite  is  the  most  stable  mineral. 

The  following  minerals  are  probably  present  in  the 
drift  based  on  x-ray  diffraction  analysis  of  drift  samples 
over  the  same  bedrock  unit  (Belly  River  Formation)  ap- 
proximately 300  km  (100  mi)  south  of  the  Sand  River 
area:  quartz,  plagioclase,  dolomite,  calcite,  very  minor 
gypsum,  montmorillonite,  illite,  kaolinite  (Wallick, 


1981).  The  results  of  the  chemical  equilibrium  com- 
putations tend  to  support  the  belief  that  chemical 
variability  of  groundwater  composition  Is  constrained 
by  equilibrium  reactions  with  the  aquifer  mineral 
assemblage.  The  carbonate  and  clay  minerals  seem 
particularly  important  in  controlling  composition 
because  the  values  of  I  obtained  were  closest  to  zero. 


Summary  and  conclusions 


The  general  hydrochemistry  of  shallow  groundwater  in 
the  Sand  River  area  is  characterized  by  TDS 
<  1000  mg/L  in  groundwater  from  aquifers  that  overlie 
the  Lea  Park  marine  shale,  while  TDS  values  generally 
exceed  1000  mg/L  in  groundwater  from  aquifers 
associated  with  the  Wapiti  Formation.  The  differences 
in  chemical  composition  between  winter  and  summer 
samplings  of  the  same  wells  plotted  on  frequency 
histograms  show  that  although  the  mean  differences 
were  statistically  indistinguishable  from  zero,  the 
distributions  of  sodium,  magnesium  and  chloride  were 
skewed  in  the  direction  of  higher  summer  than  winter 
concentrations.  This  was  attributed  to  flushing  of  these 
highly  mobile  cations  from  the  soil  and  unsaturated 
zone  during  spring  and  summer  infiltration.  When  the 
concentration  differences  of  major  chemical  consti- 
tuents were  plotted  as  a  non-scalar  bar  graph  along  the 
traverse,  positive  deflections  in  Na,  Mg,  Ca  and  SO4 
percentage  concentration  differences  tended  to  occur 
in  aquifers  on  the  western  half  underlain  by  the  Wapiti 


Formation.  Negative  deflections  tended  to  occur  in 
aquifers  in  the  eastern  half  underlain  by  the  Lea  Park 
Formation.  These  observations  were  explained  by 
noting  that  the  groundwater  flow  systems  developed  in 
the  area  where  the  Wapiti  Formation  subcrops  are 
deeper  and  have  a  longer  residence  time  so  that 
groundwater  can  become  more  mineralized.  The 
relative  effect  of  flushing  of  salts  and  groundwater  mix- 
ing on  chemical  variability  would  be  less  in  aquifers 
associated  with  the  Wapiti  Formation  compared  to 
those  underlain  by  the  Lea  Park  Formation  because  of 
the  greater  amount  of  groundwater  in  storage  in  the 
Wapiti  Formation. 

Examination  of  the  disequilibrium  indices,  computed 
by  means  of  PHCALC  and  SOLMNEQ,  showed  that 
carbonate  and  clay  minerals  were  the  most  Important 
solid  phases  entering  into  chemical  equiiibrium  with 
the  groundwater.  Chemical  equiiibrium  is  believed  to 
be  a  stabilizing  influence  on  groundwater  composition. 
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Appendix  PHCALC  Interactive  Computer 


Model  for  obtaining  equilibrium  pH, 


t3.71tjE-2,4.i6E-7,4.34E-ll,4.59£-5,-4.6312£3, 
2l.a2£3,3.55tJ.-2.9bE3, l.E-14, l.33i;4/ 
T£MPC0tA,8,CJ»**(3-298.)/C3*295..1.957)**L:3CC) 

133  WRITE(2,2J 

2  F0RM4TC1X,SP,R0GR4H  PHCALC   IS  READY   PO^  USE^.5, 

l/»  I!t»  IE»JT£R  DATA  AS   REOU 1 ED .  3  ) 

t^RITE(6,  lan 
131  rORMATClHl) 

KRITEC2,22) 
22  FORMATC/,  U,  JiNOEX  N'0.=$,1X) 

PEA0C1,2223  SAMPLE 
222  F0RHATC2A8) 

v.RITE(2,33) 
33  FORMAT(/,lX,SCALCIUM(PPM3=5,U) 

REA0C1,4)  SCA 

4  F0SMATCF13,2,' 
WRITEC2,5) 

5  FOR.IATC/,  IX,  SMAG.mESIU.H(PPM)=3,  U) 
P£A0C1,4)  SHG 

KRITEC2,6) 

6  FORMATC/,  IX,  iSODIUM(P?h')sS,  U) 
READC1,4)  SNA 

«»ITEC2.7) 

7  FORMAT  C/, IX, SPOTaSSIUHCPPmJ  iX) 
READCU-S)  SK 

WRITE(2,a3 

8  FORMATC/. lX,iCHLaRID£CPPM)=i, IX) 
REA0C1,4)  SCL 

HRI7Et2,91 

9  FORMAT  (/,  1X,SSULFATE(PFM3  IX) 
R£A0tI,4)-SSi:4 

WRITEC2,  lia) 
13  F0RHATC/,lX,lCARa0N4TE(Pp-1)=5,lX) 
READC1,4)  SC03 
-RITEC2, 11) 

11  FORMAT  C/,  lX,S!jICARl:O^ATEtPPf')s«,lXJ 
BEADC1,4)  3HCQ3 

wRITECS, 115) 
115  F0Rh4TC/,lX,JN03(PPM)=j,ix) 

R£ADCl,4)  SN03 

W&ITEC2, 116) 
115  FORMATC/, IX, SSI02(PPM) aS, IX) 

REAOd,-!)  3102 

WRITEC2, 12) 

12  FORHAT(/,lX,5FI£i.O   TEMP.    (OEG  F)  -:.U) 
REA0Cl,4)  TEMP 
XTEMP=0.5555*(TEhP-32. ) 

WRITEC2,  13) 

13  F0RM4T(/,1X,4FI£l0  PH=J,ix) 
R£ADC1,4)  PM 

KP  I  TEC  6,  3^3)  SA.SPi,£,SCA,S--G,SNi,S'<f  SC'.,SS0<i,SC0  3,SHC03#SNa3 
1,SIC2,XTEHP,PH 

333  FORM4T(/,30x, SINPUT  DATA   ECiO i , / / , 5 3 i -PlE = ■ , 2 A 6 , 2 X , $C A « 5 , 

:F7.2,2X,iHG=3,F/.2.2X,;NAij,F7.2,2<,5<:;,=-7.2,/,lx, 
1SCL=S,F7.2,2X, SS0<=5.F7.2,2X, 5C:3=3,-7.2,2x, S-C03=S, 
1F7.2,2X,5N0J=;,F7.2,/,IX,SS::2=S,'7.2.2<,?TCC):$,F7.2, 
12x,lPH(FIELJ)=i,F7.2) 
ZCA  =  SCA/4i)..-;(j£3 
Z^.G  =  SMG/24.32£3 
Z\A  =  SNi/22.5:)l£3 
Z'<  =  SK/39.  133E  J 
ZCU  =  SCL/3b.<:37E3 
2304=5304/95. 3S5E3 
7CO3  =  SC03/6i).E3 
Z-1C03  =  SHC03/61.EJ 
2'M03  =  SN03/62.E3 
ZSI02>SID2/6^.e9sE3 
BTEMP3XTEMP+273. 

BHEN  =  £XP  (TEMjJCOCHHEN,? TEMP,  AHEM) ) 

BlH  =  EIPCTEMPCOC-ilM,  F-TE'-P,  AlH)  ) 

B2HeEXPCT£MPC0(-<?i,BTEHP,  42:^)  ) 

Bw  =  EXPCTEMPCOCH.,BTEHP,Avi)) 

ESP  =  EXPCTEMPCO(.HSP.BTEMP,ASP)) 

ACTl»i1.9 

ACT2ia.7 

ZCAC03=ZMGS0iZN4S0»ZCAS0=3. 

ZCAC0=BSP/6.E-4 

ZCA=ZCA-ZCAC0 

0H=14.-PH 

Z0H=ia.*»C-0H) 

ZM»1C.««C-PHJ 

ZC03»CB2H.ZHC03)/(ACT2.ZH3 

Ki=a 

3  X5A=ZCA 

C  CALCULATE   IONIC  STRENGTH 

U«SQRTC2.«CZCAtZMG+ZS04tZC03)+2.5^CZKA*ZK+ZCL 

1+ZHC03+ZN03*Z0H+ZH+ZNAS0) ) 

Ex»-0,509«U/(l.^*l.3*U) 
C  CALCULATE  ACTIVITY  COEFFICIENT   FOR   MCSOVAL.  10^3 

ACT1=IE.*.EX 

C  CALCULATE   ACTIVITY  COEFFICIENT   FOR   DJ/AL.  lOSS 

EX.EX«4. 

ACT2«10.*.EX 
C  TREATMENT  OF  HAG.NESIU1  SULFATE   IC-N  OAIR 

lF(ZS04)7-96,  796,791 

791  IFCZMG)796,795,792 

792  AiD4CT2.*2 
Bai-C5.9E-3+AA*ZMG+AA.ZS04) 
CC'AA«ZMG«ZS04-5.9E-3«ZMG30 
XXXX. 83*33-4. a*A4.CC 
IFCXXXX)793,/93,794 

793  X3=3.a 

GO  TO  795 

794  X3=(-q9-SaRTCXXXX))/C2.B*AA) 
755  ZMGS0=ZMG50+xi 

ZMC»ZMG-X3 
ZS043ZS04-X3 


2,  and  disequilibrium  indices 


C  TREATMENT  OF  CALCIUM!  SULFATE  ION  P*IB 

79S  IFCZS04)29S,296,291 
291  IFCZCA)296,296,292 
252  A4.ACT2-»2 

BS»-(4.9E-3+4A*ZCA*AA*ZS04) 

CC«A4*ZCA.ZS04-4.9E-3«ZCAS0 

XXXXsaB.95-4.iJ«44«CC 

IF(XXXX)293.293,294 

293  X4.a. 

GO  TO  235 

294  X4.C-dB-SQRTCXXXx))/(2.a«44) 
29b  ZCAS0=ZC4SD*x4 

ZCASZCA-X4 
Z504IZS04-X4 
C  TREATMENT   OF   SGOIU**   SULFATE   10^  PAIR 

295  If CZSQ4)396, J96,i91 
301  IFCZNA)39b,395,3g2 
392  44.ACT2 


CC»»-C1.9E-l*44«ZN4*44.ZSf^4) 

CC=4a*ZN».Zsa4-l .9E-l*iN4SQ 

XXXX:ij=.3d-J.  J.44«CC 
IF(XXXX)393, J93,J94 

393  Xb=3. 

GO   TC  395 

394  X5:(-38-3QRT(XXXX))/C2.3»4A) 
393  ZNAS0=ZN430*xb 

ZNAlZMA-Xb 
Z3Q4:ZS04-X5 

C  TEST    IF   CALCIUM   CO'KENTR  A  T I  ON    IS  CONSTA.ST 

396  K1=K1+1 

ERR0Ri(XCA-ZCA)/XC4 

IFCABSCERROR J  .GT. .21 )   GO   TO  3 

ZH«ZH/4CTl 

EOH=  C82H/aSP)*ZCA*ZHC03*ACT2 
EOPMa-ALOGiacEOH.ACTl) 

SIC>AL0G13CCACT2.ZCA«ZHCO3*B2"»ie.««P-)/35') 

SI0>AL0Gl3(CACT2.4CT2«4CTl*ACTl»ZCA.Z-^G»ZHC03.Z-(CO3 

l«B2wB2H»ia.««(2.«PM))/2,34E-l7) 

SIG:AL0G13CZCA*Zsa4*ACT2»ACT2)*4.61-2.49E-3«XTE-P+ 

15.92£-5»«TEmP*XTEMP 

ZC03=SS?/(ZCA«CACT2«4CT2)) 

ZHC03:i(ZH.4Cr2«ZC03)/B2H 

Zm2CG3=(ZM*ZHC03.(4CT1«4CT1))/31H 

PC02»Zh2C03/BM£n 

CSP=ZCA.ZC03«(ACT2*4CT2) 

C1H={ZH.2HC03*(ACT1*ACTI) )/CZH2CC3) 

C2H=CZH«ZC0J»ACTl)/(ZHC03) 

TCA=ZCA+ZCA30 

STC023ZH2C03*ZHC03tZC03 

4IABK34L0Gia(ZNA*4CTl)+2.»AL0Gl3(ZSIC2)*=-i-C . 34 6* . 33323 • TEMP ) 
AMKiALOGlJCZ**  AC  Tl)->'2.«4L0Gia(ZSI02)t'H-C.;i529»rEHp. 2.931) 
AA.MKiAL0Gie(ZCA*4CTn+PH-ci8.84-.371lJ.TE>i9) 

AIK  =  .5*4LQG1  -J(ZK.4CT1)*.25«AL0GI£CZ'G«4CT2)  +  1 .2«AL0GIB(ZSI02) 
1        +1 . l*PH-( 1.692-.2a5a56«T£MB) 

ACAfK  =  .l666.(4L0Gli«CZCA.4CT2)*S.«ALjG:aCZSI02)*2.*OH-C.3U43* 
1  TEMP-16.6)) 

AKGr2.»ALOGl<3(ZSlO2)-(.02*TEMP-8.75) 
C  CONVERT  MOLALITY   TO  PPh  » 

ZCA  =  ZCA.4/;.  Ja£3 
ZMG  =  Z>1G»24.32£3 
ZNA=ZNA.22.991E3 
ZK  =  Z'<*39.13CE3 
ZM2C03=ZH2C03'62.£3 
ZC03=ZC03«63.£3 

zcAcoszCACo*i3a,a8E3 

ZS04sZS04«96.366£3 
ZHCG3=ZHC03»61 .£3 
ZMGS0=ZMGS0»1 10.386E3 
ZNASO  =  ZNA3O*l20.;;56E3 
ZCA3C=ZCAS0«l36.14o£3 
ZCL»ZCL*35.457E3 
fc-RITE(6,14)  ZCA,ZMG,ZNA,ZK 
HRITE(e, 15)  ZH2C03,ZHCC3,ZC03,ZSC<.ZC. 
WRIT£C6,15)    ZCACO,Z*<GSO,ZNASO,  ZCAS3 
WRIT£C6,17)    EQPh,    ACTl,    ACT2,    SIC,  Si: 
».RITEC6,232)  PC02,SIG 
2-^2  FGRMir(/,lX,iPCO2  =  S,fcl0.2,2X,SSiTG'"i£,£13.2) 

14  FORy.ATC/,  lX,iCALCULAT£C  CATIOv  C  C  •■<C  £^  "3  i  T I C  "•.  3  (  b  =  h.  )  s  , 
:/,lX,JCA=5,FlC,2.2X,iMG=S,F13.2,2X,5Si:5,ri:,2,2X,$K»S,Fie.2) 

15  FOP  1ATC/,1X,SCALCULATE0   CARBOSiTE   iNj  asICa  CC^CEN TRA riONS 
l(PPM)S,/,lx,iH2C03-S,Fl3.2,2X,S-CCJi5,FlJ,2,2x,5Cnj:5, 
2FlP.2,2X,JS04=S,Fia.2,/,lx,5CL»5.=ir.2) 

1*  r.-RH4TC/,  IX.iCALCJLATEC    ION  Pil=  CC  ^Ci'- r  ;  i  t  i  c-,s  ( ?pu )  s  . 

W.  IX,  ?C4C03  =  5,F17.2,2X,  3-GS04rS,F:  J,2,2X,  SMS:i  =  S,F13.2, 

1?., JCiS0-=S,F1^.2) 
17  F0=  1AT(/, IX, jEQUIL  PHsS,F5.2,3X,Si:"i:3,-5.2,3s,iACT2« 

13,F5.2,3X,S3lCsS,El?.2,2x,JS:D=S,i.;r.2) 

*-rtIT£C6,23)  4lABx,AHK.4AN<,AIK,iCA-<,-<C 
23  FORMiTt/,  lX,iICAa-K):(   ;  ,  E  1  3  .  2  ,  I  X  ,  5 1  C  K --K  )  =  3  ,  £  1  :  .  2  ,  1  X  ,  S  IC  AN-K  )  s 

1   5,Eie.2,//,lX,SICIL-K)=r$,E12.2,:x,3ICC>'--<>a   i  ,  £  1  2  .  2  ,  1  X  ,  S 1 C  K- 

1G)=  i,£13.2; 

»^IT£C2, 19) 
IS  FOR«ATt/,lX.  JANY   J10RE  0ATA7S) 

READC1,19)  ASS 
19  F0RmaT(A3) 

IF(4NS.E0.3m»'ES)    30   TD  100 
999  STOP 

fcsD 


N.L.C.  -  B.N.C. 
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